Observation of Spin Glass Dynamics in Dilute LiHo ;c Yi_ x F4 
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AC susceptibility measurements are presented on the dilute, dipolar coupled, Ising magnet 
LiHo a; Yi_ I F4 for a concentration x = 0.045. The frequency and temperature dependences of the 
susceptibility show characteristic glassy relaxation. The absorption spectrum is found to broaden 
with decreasing temperature suggesting that the material is behaving as a spin glass and not as an 
exotic spin liquid as was previously observed. A dynamical scaling analysis suggests a spin glass 
transition temperature of 42.6 mK ±2 mK with an exponent zv — 7.8 ± 0.23. 



Spin glass behavior is an effect resulting from quenched 
disorder that has been extensively studied for many years 
and is largely well understood. It was therefore surprising 
that the system LiHOj.Yi_3.F4, & seemingly ideal dilute, 
dipolar coupled Ising magnet, would exhibit an unusual 
spin liquid state at a low concentration of magnetic ions, 
x = 0.045 [H, Q. This spin liquid, or "antiglass" phase 
was most notably characterized by a narrowing of the 
absorption spectrum }("{lo) as temperature was lowered. 
This was found in stark contrast with a higher concen- 
tration (as = 0.167) spin glass state where the absorption 
spectrum broadens [3, 0, Q and with theoretical predic- 
tions that the spin glass state should persist to concen- 
tration in such a system with long-range interactions [H . 

Later publications found numerous interesting effects 
in the spin liquid state including very narrow absorp- 
tion spectra (narrower than what can be accounted for 
with models of glassy relaxation) with strong asymmetry, 
ringing magnetization oscillations, sharp features in the 
specific heat and a T -075 power law in the dc limit of 
the susceptibility 0, 0] ■ 

Recently there has been a large amount of activity on 
this series of materials, both theoretical and experimen- 
tal. In particular there is an ongoing debate on the ex- 
istence of a spin glass transition at any value of x in the 
LiHo x Yi_ l; F4 series @, _], [13, [ll[ . Much research has also 
taken place attempting to understand the effects of trans- 
verse magnetic field on the purported spin glass state (see 
for example 0, _3, [H, [_1 ) . 

In this Letter we will focus on the low-concentration 
limit of these materials (x — 0.045) in zero field. Wc 
have measured the ac magnetic susceptibility of this sto- 
ichiomctry in the hope of reproducing the exotic physics 
that was observed previously. We will show that, in fact, 
the susceptibility of this material behaves much more like 
that of a spin glass and that we are unable to reproduce 
the unusual antiglass phenomenology. 

For this work, a magnetometer based on a dc super- 
conducting quantum interference device (SQUID) was 
developed, chosen for its sensitivity and flat frequency 
response at very low frequencies. The magnetometer con- 
sists of a 375-turn NbTi primary coil wrapped on a phc- 
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FIG. 1: AC susceptibility of the x = 0.045 sample showing in- 
phase x'(f) an d out-of-phase x"(/) components. The spectra 
were obtained at temperatures 77, 80, 83, 85, 87, 90, 95, 100, 
110, 120, 140, 150, 170, 200, 230, 260, 300 and 350 mK from 
left (blue) to right (red). 



nolic form surrounding a niobium, 2nd-order gradiomc- 
tcr. The secondary forms one branch of a continuous 
superconducting loop or flux transformer which couples 
flux to a dc SQUID [2]. The flux transformer makes 
use of the Meissner effect meaning that our signal is not 
proportional to the frequency of excitation as it is in a 
standard inductive susceptometer. This is crucial for at- 
taining the low frequencies necessary to study such glassy 
dynamics. An additional trim coil is used in parallel with 
the primary to remove misbalance of the gradiometer. 

The SQUID is contained within a superconducting lead 
shield and is run in feedback mode. The entire magne- 
tometer apparatus is contained within a superconducting 
lead shield and the cryostat itself is surrounded by a lead- 
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coated radiation shield and two /i- metal shields. 

Measurements were performed on a sample of 
LiHo 2; Yi_ 2 ;F4 with x = 0.045. The sample was glued to a 
sapphire rod, the other end of which was heat sunk to the 
mixing chamber of a S.H.E. dilution refrigerator. Previ- 
ously, the specific heat of this same sample was measured 
and the sample was characterized with various X-ray 
scattering techniques as discussed in Ref. [HI . For most 
measurements, the sample was cut to be needle-shaped 
(dimensions 0.57 mm x 0.77 mm x 7.7 mm) to reduce de- 
magnetization effects with the long dimension along the 
c-axis (also the direction of the applied field) . The result- 
ing demagnetization factor is AttN = 0.493 [T3|. Some 
measurements were reproduced using the same sample 
cut to a different aspect ratio (0.57 mm x 0.77 mm x 3.31 
mm) with demagnetization factor 4irN = 1.104. Match- 
ing the results for two differently shaped samples allowed 
us to calibrate our magnetometer and confidently deter- 
mine the correct demagnetization correction. 

Both detailed frequency scans at constant temperature 
and detailed temperature scans at constant frequency 
were obtained. It was carefully checked that the sam- 
ple was in equilibrium for all measurements by waiting 
several hours or more before measurement at a given tem- 
perature and also by taking multiple spectra to check for 
reproducibility. The applied ac magnetic field was kept 
below 20 mOe at all times to ensure that there was no ap- 
preciable heating of the sample. A range of fields around 
this value was tested and none were found to cause any 
significant difference in the spectra that might indicate 
a heating effect. Measurements in the frequency range 
0.001 Hz to 2 kHz are presented here. The accessible 
frequency window of these experiments was limited at 
the high end by frequency dependent background signals 
and phase shifts, and the SQUID feedback electronics. 
At the low end, it is fairly impractical to perform any 
measurements appreciably below 1 mHz. 

The measured frequency scans of the complex ac sus- 
ceptibility are typical of glassy relaxation, showing broad 
peaks in x"^) an( i suppression of x'{<jj) at higher fre- 
quencies as shown in Fig. [I] The characteristic time con- 
stants of relaxation, which may be parametrized by the 
peak frequencies in \" , /Max, are shifted to lower fre- 
quency with decreasing temperature as thermal fluctua- 
tions become weaker and are unable to excite the system 
over energy barriers. 

Several studies of the x( w ) spectra of spin glasses 
have been pcrformed_prcviously including work on the 



and on other materials 18, [1 



LiHozYi-a^ serie 

l20l |. Despite being quite broad, the absorption spectra 
measured here are in fact narrow when compared to some 
other spin glass measurements [1, 18|. Similar glassy re- 
laxation has also been observed in t he g eometrically frus- 
trated spin ice material Dy 2 Ti 2 07 [2lL l22T|. 

Normalizing x"( w ) by the peak height x'^iax ana ^ di- 
viding the frequency by the peak frequency /Max super- 
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FIG. 2: Frequency scans of \" f° r various temperatures nor- 
malized by the peak height XMax on the vertical axis and 
by the peak frequency /Max on the horizontal axis. A clear 
broadening of the absorption spectra can be seen as the tem- 
perature is lowered from 260 mK to 77 mK. The full width at 
half maximum (FWHM) in decades is shown in the inset as 
a function of inverse temperature. The FWHM can only be 
obtained at a limited set of temperatures where the experi- 
mental frequency window envelops both sides of the curve. 



imposes the absorption spectra. The result is a clear 
broadening of the spectra with decreasing temperature 
as seen in Fig. [2l The full widths at half the maximum 
of the spectra (FWHM) are plotted in the inset of Fig. [5] 
The FWHM appears to level off at ~ 1.4 decades at 
higher temperatures. Though we have not shown fre- 
quency scans here, an 8% sample also showed broaden- 
ing of the absorption spectra with decreasing tempera- 
ture. This behavior is qualitatively similar with that of 
a higher concentration, x — 0.167, material in this se- 
ries [2j. 

At the higher temperatures studied, the absorption 
spectra show low and high frequency limiting behavior of 
x" ~ ljJ 1 and x" ~ o->~ ' 75 respectively. This is not con- 
sistent with the Debye model with a single time constant 
where the limits are uj 1 and uj~ 1 . As the temperature is 
reduced and the curves become broader, the tails of the 
spectrum become less steep. At the lowest temperatures 
studied here, the high frequency limit can be seen to be as 
shallow as x" ~ w~ 63 . These power laws constrain what 
functions might fit the data. Several standard fitting 
functions including the Debye model with various distri- 
butions of relaxation times, the Davidson-Cole form [23| 
and the Havriliak-Negami form [24[ were considered but 
were not able to fit the full frequency range studied. 

An important aspect of glassy systems is the critical 
slowing of dynamics or a divergence of the characteristic 
time constant of the system r. The most robust way of 
determining r would be to take the limit of x" /^x' as 
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lu — > [25[, but such a limit is only achievable in a very 
small temperature window due to the broadness of the 
spectra and experimental time limitations. Thus, in this 
work, we have chosen to parametrize the dynamics of the 
system with T Max = l/27r/ M ax- 

At first glance, these peak positions appear to roughly 
follow an Arrhenius law, Tuax(T) = tqa exp (— Ea/UbT), 
at least at higher temperatures, as shown in Fig. [3ja), 
giving Ea = 1-57 K and t a = 0.32 fj,s. However there 
is some curvature thus the Arrhenius fit is only able to 
accommodate the higher temperature data points and 
significant deviation occurs at lower T. A likely sce- 
nario at lower T is a dynamical scaling law of the form 
TMax = Tn(T/T„ — X)~ zv , that is predicted to apply to 
spin glasses [25|. Such a fit, plotted in Fig.[3l[b), is quite 
successful for the data below ~ 200 mK giving a transi- 
tion temperature T g = 42.6 mK ±2 mK and an exponent 
of zv = 7.8 ± 0.23, very near the exponent 7.9 obtained 
in Monte Carlo simulations (2^]. Above 200 mK, as we 
move further from T g and out of the critical regime, the 
power law behavior breaks down and appears to give way 
to an Arrhenius law. 

The overall intrinsic time scale of this material, tq ~ 
20 s, extracted from the power law fit, is extremely long 
(in Euo.4Sr . 6 S, for example, r ~ 2 x 10~ 7 s [26J) and 
explains why obtaining equilibrium data anywhere close 
to the transition temperature becomes completely im- 
practical. For example, measurements with a frequency 
of ~ 3 x I0~ 5 Hz would be required to properly study 
even T = l.5T g . This difficulty would also apply to other 
measurements such as nonlinear susceptibility. The near- 
Arrhcnius law that is observed is likely a consequence of 
measuring at temperatures far above T g . 

Detailed temperature scans at constant frequencies 
were also obtained and x'(D 1S shown in Fig. [4j Pre- 
dictably, one sees an increase in \' with lower tempera- 
ture until the frequency of relaxation of the sample be- 
comes slower than the probe frequency. At higher tem- 
peratures, where the sample's time constant is fast, the 
susceptibility measurements shown can be considered to 
be in the limit of static susceptibility. However, below 
approximately 90 mK, the sample's relaxation is so slow 
that even 1 mHz cannot be considered to be in the dc 
limit and there is a downturn in x' ■ Nevertheless above 
100 mK or so, we can effectively probe the static suscep- 
tibility of the material. 

Comparison of these results with those obtained by Re- 
ich et al. [H, 0] shows a fair bit of agreement in certain 
respects. In both cases the peak frequencies roughly fol- 
low an Arrhenius law and match quantitatively as shown 
in Fig. [3ja) . The static limit of the magnetic suscepti- 
bility shows somewhat similar behavior as a function of 
temperature (see Fig.|4j. However, our results show that 
the widths of the absorption spectra do not narrow but 
rather broaden with lower temperature, a result more 
consistent with a spin glass. In other words, we do not 
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FIG. 3: Different scenarios for scaling of the time constants 
obtained from the peak frequencies of the absorption spectra, 
i"Max(T). (a) Data from this work fitted by an Arrhenius law 
at high T (straight line). Also shown are the time constants 
TMax from Ref. [IJ (dots) and Ref. [(J (crosses), (b) Data from 
this work (circles) fitted by a dynamical scaling law at low T. 
(c) Residuals and estimated error bars for the Arrhenius and 
scaling law fits. 



see "antiglass" physics reported in Ref.'s [1,0]. 

Papers published more recently 0, 0] show results that 
differ strongly both from work by the same research 
group [3, 0| and from our results presented here. The 
shape and width of the absorption spectra are qualita- 
tively different in Ref. [6j from the results in this work. 
The peak frequencies also do not match with Ref.'s [H, 0] 
and clearly do not follow an Arrhenius law, as shown 
in Fig[3ja). Additionally, Ref. 0] states that the static 
susceptibility obeys a power law T~ a with an unusual 
a = 0.75, a much shallower temperature dependence 
from that seen in this work and by Reich et al. [H, 0] (se e 
Fig. [4]) . Previous work by the authors of this paper [l6( , 
measuring the specific heat of several stoichiometries in 
this series, found consistent, smooth features also in dis- 
agreement with Ref.'s where unusual, sharp features 
were observed. 

The ac susceptibility of L1H00.045Y0.955F4 was also 
measured recently by a third research group Q and a 
third, distinct temperature dependence is observed. At 
higher temperatures, where the response is fast, Ref. [§| 
matches well with our results here. At lower tempera- 
tures, the discrepancy is likely a result of the magnetic 
field sweep rates used in Ref. corresponding to fre- 
quencies of measurement that are too fast to be con- 
sidered in the static limit and leading to an apparent 
reduction of \ ( scc Fig.@|. 
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We also compare our results to very recent Monte Carlo 
simulations of this material [27j and find that there is in- 
deed very close agreement with the static limit of the sus- 
ceptibility as a function of temperature (Fig. QJ. There 
is also qualitative ag reement between the specific heat 
measured in Ref. [l(| and that calculated in the simu- 
lations [27j. These and other recent Monte Carlo sim- 
ulations [1, 27, find no divergence of the spin glass 
susceptibility suggesting that there is no finite tempera- 
ture spin glass transition even if this system is viewed as 
a perfect Ising model. 

However, classically, the mean-field theory of Aharony 
and Stephen [f| would lead us to believe that there should 
be a spin glass transition all the way down to x — 0. It 
has been suggested by Ghosh et al. 0] that off-diagonal 
terms inherent in the dipolar interaction can introduce 
quantum fluctuations, leading to a spin liquid state. 
Recent theoretical work 14j analyzing a more detailed 
model incorporating random fields resulting from the 
hyperfine interaction and off-diagonal dipolar coupling- 
maintains that quantum effects are not strong enough to 
stabilize a spin liquid state and that the system should 
therefore undergo a spin glass transition. They estimate 
a T g of roughly 35 mK, quite close to the 42.6 mK that 
we have determined in this work. Additionally, Schechter 
and Stamp [3] predict a significant slowing of the dy- 
namics as the concentration x is lowered to 4.5% as a 
result of the nuclear hyperfine coupling. 

In conclusion, our measurements of LiHO; C Yi_. 1 .F4 have 
not shown the exotic antiglass physics that was observed 
previously 0, 0, 0, 0|- Instead, the absorption spec- 
trum broadens with lower temperature, consistent with 
behavior expected of a spin glass. A dynamical scal- 
ing analysis using t Max(^) obtained from the maxima 
in x" implies the existence of a spin glass transition 
around T g = 43 mK. While the temperature depen- 
dence of the static susceptibility (Fig. and specific 
heat [lj| of this sample are similar to the Monte Carlo 
simulations [27|], dynamical scaling provides compelling 
evidence for a finite temperature spin glass transition, a 
conclusion which is supported by the theoretical work of 
Ref. 14| . An extremely slow response of the system has 
also been observed and illustrates clearly the need for 
very low frequency measurements when studying dilute 
LiHo x Yi -xF4 and other similar systems. The ongoing 
debate (H, 0, 0, EH on the existence of a spin glass tran- 
sition at higher concentration LiHoo.i67Yo.833F4 might 
also be resolved through more careful attention to the 
diverging time scales involved. 

We have benefited greatly from discussions with 
M. J. P. Gingras. Thanks also to S. W. Kycia and 
A. Gomez for discussions and help with sample charac- 
terization. Funding for this research was provided by 
NSERC, CFI, MMO and Research Corporation grants. 



o 
O 



H 0.1 

u 



8 



0.001 Hz 
0.01 Hz 
0.1 Hz 
1 Hz 

from Ref. [9] 

from Ref. [2] 

from Ref. [27] 
y-0.75 



0.1 



Temperature (K) 



FIG. 4: Temperature scans of \' a t four different frequencies 
of excitation. We show for comparison experimental and the- 
oretical results from other research groups: the DC limit of x 
from Ref. [3 (as published) , x' taken with a linear field sweep 
from Ref. [SJ| (after correcting for demagnetization), Monte 
Carlo simulations [27[ (arbitrary units, scaled to match at 
high T) and a T -0,75 power law that was proposed in Ref. 0. 
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